1 



QCD in Nuclear Collisions* t 
A.H. Mueller a 

a Department of Physics, Columbia University 
New York, New York 10027, USA 

The relationship between RHIC and HERA data is explored using the idea of saturation 
(color glass condensate) as a unifying framework for interpretation. A description of the 
early stages of a heavy ion collision is given consistent with the RHIC results. The relevant 
QCD dynamics at the various stages of a heavy ion collision, from production of gluons 
to equilibration, is reviewed. Some comments on current phenomenology are given. 

1. Introduction 

This lecture is devoted to a general semiquantitative description of the early stages of 
a heavy ion collision. This is a rapidly evolving subject and it is not yet clear how good 
a description can be given in terms of QCD. 

In the earliest stages of the collision gluon occupation numbers are estimated to be 
moderately large suggesting a classical gauge theory description. This fits in nicely with 
the saturation picture (color glass condensate) of a high-energy heavy ion wavefunction. 
The parameters necessary for such a description seem compatible with an analogous de- 
scription of small-x and moderate Q 2 HERA data. 

Various discussions of RHIC data in terms of gluon saturation, which dicussions have 
much in common, are reviewed. The general picture appears to be coherent although per- 
haps not yet compelling. The reader should be warned that many of the numbers quoted 
in the text are theoretical estimates. There are large uncertainties present, neverthelelss, 
I believe the estimates to be reasonable. 

2. General features of the early stages of high energy heavy ion collisions 

2.1. Gluon number densities at the initial time 

Consider a centrai collision at RHIC, say at ^/s = 130GeV. The measured transverse 
energy is ~ 500GeU|l[| at centrai rapidities. Since some of the initial transverse 
energy goes into longitudinal flow one certainly has 



a T > 500GeV 



dy 
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for the initial transverse energy distribution, which we assume is predominately in gluons. 
Thus one can write 



dE™ 
dy 



dN l „ n 
9 -AK 



dyd 3 b 



A g ■ K > ÒOOGeV (1) 



with b±, and b z the coordinates of a point in the overlaping nuclei and where Ab z is the 
longitudinal width of the volume occupied by gluons at their time of production, A g is 
the area, transverse to the colliding beams, occupied by the gluons and K is a typical 
gluon transverse momentum. Taking Ab z — 1/K and A g ?» 150/m 2 gives 

dN in 

and with < K >^ 2/3 (see later) 
dNÌ n 

where (2) and (3), as well as < K 2/3GeV represent average values, averaged over 
ali transverse coordinates. At the center of the nucleus, bj_ = 0, the more appropriate 
value of K is about 3/2 the average value or K ~ Q s (b± — 0) ~ IGeV, with Q s the gluon 
saturation momentum, so that 

Equations (2-4) give a reasonable guess as to the initial gluon densities averaged over the 
production region (Equations (2,3)) and at the center of the production region (Equation 
(4))- 

2.2. Occupation numbers 

The initial gluon number densities discussed just above can be converted to gluon 
occupation numbers (phase space densities) by 

(2tt) 3 dNì n (2tt) 3 dNì n 

Jg 2 • (N£ - 1) d 3 pd 3 b ~ 2(N 2 C - 1) dyd 2 b ± d 2 p ± y } 

where we use Ay = Ap z /p z ~ Ab z Ap z and where the 2 in the denominator refers to the 
number of spin states available to gluons. We further make the approximation 

dN T „ 1 dN ? (6) 



dyd 2 b±d 2 p± irQ 2 dyd 2 bj 
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Using (4) and (6) in (5) gives 

ff ~ 1.5 (7) 

at 6j_ = 0. Thus, it would seem that gluon occupation numbers are moderately large at 
initial times thus fitting well the picture of strong gluon field dominance at early times in 
heavy ion collisions. 

2.3. The nuclear wavefunction 

It is convenient to view a head-on (centrai) ion-ion collision in a frame where one of 
the ions is at rest and the other carries a large longitudinal momentum. We view the 
scattering in a light-cone gauge which furnishes the smoothiest connection between gluon 
quanta present in the wavefunction and the gluons produced at the time of the collision. 
Roughly, produced gluons at mid-rapidity are already in the wavefunction of the fast ion 
and are merely "freed" by the collision |2|||. The picture is illustrated in Fig.l. The initial 
wavefunction has a large occupation number for gluons, f g ~ 1/at, which results in a high 
density of produced gluons at the time of production. The fact that f g is large in the 
high-energy ion's wavefunction is one of the features that led McLerran and collaborators 
to cali it a color glass condensate. 




Figure 1. 



2.4. The McLerran- Venugopalan model 

The McLerran- Venugopalan |Q] model is a simple but maybe not unrealistic model of the 
small-x gluonic part of the wavefunction of a heavy ion at RHIC energies. The Feynman 
graphs contributing to the wavefunction|| are illustrated in Fig.2. The gluons in the 
wavefunction correpond to a quasiclassical picture because only tree graphs are involved 
even when one squares the wavefunction to get the gluon number density. The gluon 
number density takes a very simple forni 
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Figure 2. 



dN g dxG a iV c -l f d 2 x± 



dyd 2 bd 2 k ± d 2 bd 2 k ± An A aN c J x 
where the saturation momentun, Q s , is given by 



r f^ e -i k± . x ^ 1 _ eX iQV4 ) (8) 



Aix 2 aN c 



x ) = 1 ^—Y^v^-^p^G nucleon (9) 

as a function of the impact parameter in the nucleus and of the Bjorken variable, x. 
Calling 



^ = 2VW^¥pxG nucleon (10) 
one can write the saturation momentum as 



2 _ AnaN c dxG a nn 

^ s ~ ìv c 2 - 1 d 2 b ■ 1 J 

It is straightforward to get the gluon occupation number in the wavefunction from (8), 
using a formula analogous to (5), and one finds 



f 9 = 4r f ^e~ ik —{l - e-*ì^) ~ 4r r -e- tfe -^ (12) 
' 9 aiV c y ttx 2 ± v ' aN c J 1 t y ' 

where the approximate evaluation, given by the t-integral, is valid when k 2 ± /Q 2 is not too 
large. This form of the integrai neglects a weak x\— singularity which is present in Q 2 
via the scale of xG in (9). 

From (12) it is apparent that f g ~ 1/a when k 2 ± /Q s is not large. Occupation numbers of 
size 1/a are characteristic of classical field theory solutions, here using the valence quarks 
of the nucleons of the nucleus as the sources of color charge. The McLerran-Venugopalan 
model gives the gluons in the wavefunction of the nucleus as the Weizsàcker- Williams 
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fìeld generateci by the valence quarks of the nucleus. At RHIC energies this may be a not 
bad approximation. At LHC energies quantun evolution in the longitudinal momentum 
(BFKL evolution) will also be important thus further increasing the gluon numbers by 
raising the saturation momentum, Q s |||9|]. For gluons having kj_/Q s on the order of 
one, which gluons will dominate the initially produced energy in an ion-ion collision, the 
occupation number f g ~ 1/oc. The fact that f g never grows larger than a fixed Constant 
times [aiV c ] _1 when k±/Q s ~ 1 is the essence of the idea of gluon saturation [|T0| , pA| and 



corresponds to the dominance of the initially produced energy density by fields F^ u ~ 1/g. 
3. Relationship with HERA phenomenology 



3.1. The Golec-Biernat Wiisthoff model |l2fl 

Many of the properties of the small-x wavefunction of a high-energy heavy ion are also 
present in high-energy protons and have been extensively studied at HERA. There is some 
evidence for saturation at HERA, although most of the evidence is indirect and makes 



use of models of the Golec-Biernat Wiisthoff type ||T2|--[16 1 . 

In the Golec-Biernat Wiisthoff model one views virtual photon-proton scattering in 
terms a quark-antiquark dipole, coming from the 7*, scattering on the proton. QCD 
evolution is put into the proton so that the scattering is that of a bare dipole, having no 
evolution, colliding with a highly evolved proton wavefunction. The situation is illustrated 
in Fig.3 while the formulas are]T2| 




Proton 



Figure 3. 



Q 2 



/d 2 x± / dz 
Jo p 



(13) 



and 
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F ? = T^~ [ d 2 x ± f dzY,^F\M^^^Q)\ 2 \a Q {l - e-^/'f + qqg part 

Oi ern J Jq Z 

(14) 

for F 2 and the diffractive part of F 2) respectively. 

The various factors in (13) have ready interpretations. The — is the factor con- 
necting F 2 to the scattering cross section of an (off-shell) transversely polarized virtual 
photon on a proton. ip T is the lowest order wavefunction for a photon of virtuality Q to 
go into a quark- antiquark pair having transverse coordinate separation x± and with longi- 
tudinal momentum fractions z and 1 — z for the two members of the dipole. e~ x ^ Q ^ 4 = S 
represents the impact parameter averaged S-matrix for the elastic scattering of a dipole 
on a proton while o"o represents the size (area) of the proton. 

There are three parameters in the model 

o"o = 2.3mò, xo = 3xl0~ 4 and A = 0.3 (15) 

where 

Ql = (xo/x)\ (16) 

in units of GeV 2 , with Q 2 S having the interpretation of an impact parameter averaged quark 
saturation momentum of the proton where the gluon and quark saturation momenta are 
related by 

Q 2 s = ^Q 2 s = WÌ (17) 

The interpretation of (14) goes much the same way There is also a term involving the 
scattering of a quark-antiquark-gluon system on the proton, but this introduces no new 
parameters. 

With the three parameters of the model given as in (15) the G-B W model does a good 
job in fitting the moderate Q 2 and low x HERA data. Saturation plays an important role 
because the dipole-proton S-matrix 

S = e x2 ^/ 4 (18) 

becomes zero when x± > 2/Q s thus limiting the contribution of large dipole configurations 
in (13) and (14) to that allowed by unitarity. This limitation is especially important in 
diffractive scattering which in perturbation theory has a strong divergence coming from 
large dipole sizes. 

In the HERA range the value of Q 2 from (16) is between 1 and 2 GeV 2 giving Ql ~ 
(2 — A)GeV 2 . As we shall see later this value may be a little large and a more realistic 
model can be expected to give somewhat lower values of Q 2 . 
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Figure 4. 



3.2. Going to impact parameter space 

One can get information on the S-matrix for dipole-proton scattering as a function of im- 
pact parameter by using 7* proton — ► J/ip proton and 7* proton — > p proton quasielastic 
scattering data]l7|. The process is pictured in Fig.4 for the production of longitudinally 
polarized p-mesons at momentum transfer A. To analyze the process one also needs a 
model of the p-meson wavefunction which is taken from earlier phenomenological anal- 
yses|TE|-|2"0|]. In the approximation that the dipole-proton scattering amplitude is purely 
imaginary one can write the amplitude in terms of the square root of the differential cross 
section. It is straightforward to arrive at 




S(b, x ± ,x)=l- / cPAe- iA \/ ^ (19) 



where 



N = (rl> PL ,il>r) (20) 

is the overlap of the wavefunctions for a longitudinal virtual photon with a longitudinal 
p— meson. S depends on the impact parameter, b, on the Bjorken x— variable, and on the 
transverse size of the dipole pair, x±, which is fixed, approximately, in terms of Q. 

In order to interpret (19) in terms quantities we have used in the Golec-Biernat Wùsthoff 
model write S as in (18)||17jl. (This form for S is not generally true, but it is not unrea- 
sonable when gluon number densities are large.) From fitting the data at x — 5xl0~ 4 and 
Q 2 = 3.5GeV 2 (x± ~ 0.4/m) one gets 



Q 2 s (b = 0) ~ 0.5Geì/ 2 or Q 2 s (b = 0) ~ l.lGeV 2 . (21) 

Although the values given in (21) have considerable uncertainity, mainly from deciding 
what value of x± to assolcate with a given Q 2 , I think they are reasonable. They are 
somewhat smaller than those that come from the Golec-Biernat Wùsthoff model. Using 
(11) to write 
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Q 2 s (b,x) 
we are led to 



2-n 2 a dxG 

N r In 



(22) 



-^-|&=o ^ 5/fm 2 and xG ~ 6 (23) 

for the values of x and Q 2 quoted above. 

3.3. From HERA to RHIC 

Now let's be really adventurous and scale (23) back to RHIC energies, x ~ 0.01, by 
using the scaling law 



Ql(x) = (^f 3 Q^o) (24) 

x 

found in the Golec-Biernat Wùsthoff analysis. This gives 



xG\rhic ^ 2.5 (25) 

a little larger than most perturbative QCD fits, but not a crazy number at ali. 

Now that we have xG for the proton at RHIC energies we can use (9) and (10) to get 
the corresponding gluon density and saturation momentum. For Au at RHIC energies 
one finds from (9), (10) and (25) 



Q 2 s \ b=0 ~ 1.25GW 2 and -— % =0 ~ 5/fm 2 (26) 

ero 

not too different from our estimate given in Sec.2.1 for the values of the typical momentum 
and number density of gluons produced in a heavy ion collision. We shall come back a 
little later to the question of the relationship between the number density of gluons in the 
wavefunction of an ion just before collision and the number density of initially produced 
gluons. 

4. What dynamics can (should) one use? 

For the wavefunction of a heavy ion, and at RHIC energies, the McLerran-Venugopalan 
model should be adequate. However, by LHC energies, gluon evolution will certainly 
become important0,[| and may even be significant at RHIC energies. The density of 
gluons in the wavefunction (26) and initially produced (2-4) is so large at RHIC that the 
idea of viewing the wavefunction and the production process in terms of hadrons and in 
terms of nucleon-nucleon collisions, respectively, does not make sense. The collision must 
be viewed in terms of quarks and gluons. 
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The McLerran-Venugopalan model uses classical gluon dynamics with random color 
sources given by the valence quarks. Because gluon occupation numbers are reasonably 
large the freeing (production) of gluons in an ion-ion collision should also be amenable to 
a description in terms of classical gluon dynamics. So far there have been two important 
approaches to this problem. 

Krasnitz and Venugopalan |2"T| , |2"2"| | have developed a numerical procedure for calculat- 
ing the early stages of a heavy ion collision in terms of discretized classical Yang-Mills 
equations. The results of these calculations will be discussed shortly. Kovchegovp3| has 



tried to calculate the gluon production cross section by looking at Feynman graphs. His 
procedure for gluon production in proton-nucleus and nucleus-proton collisions gives re- 
liable results. (Nucleus-proton scattering is not manifestly the same as proton-nucleus 
scattering in Kovchegov's calculation because of the use of a particular light-cone gauge.) 
It is stili an open question whether his procedure is valid for nucleus-nucleus collisions, 
although the results which come out are quite reasonable. 



4.1. Krasnitz, Nara, Venugopalan calculation ||22| 

Starting with a McLerran-Venugopalan wavefunction as an initial condition Krasnitz, 
Nara and Venugopalan numerically solve discretize Yang-Mills equations as the nuclei in 
a heavy ion collision overlap. Gluons appear to be freed on a time scale r ~ 3/Q s . Just 
before the collision the distribution of gluons in one of the nuclear wavefunctions is given 
by (8) and with gluon occupation numbers given by (12). Just after the gluons are freed 
the distribution of liberated gluons fits well to the empirical formula 



dN T _n 2 c -i i o.ii (27) 



dyd 2 b ± d 2 p± Att 3 aN c e ^pl+m 2 /T eff _ ì 
with T e ff ~ l.lQs and m ~ 0.08<5 S giving 

f in ~ — f28ì 

Jg ~ aN c eP^ff - 1 K ' 

so long as p±/T e ff is not too small. 

Comparing (8) to (27) one sees that c, the fraction of gluons freed in the collision, 
obeys c ~ 1/2. In the wavefunction < p±/Q s >~ 0.6 while in the initial distribution 
(28), < p±/Q s >— 1.5 so that the freeing process seems to create a lot of additional 
transverse energy which is not present in the ion wavefunction. One worrisome result 
of the calculation is that the occupation number, given by (28), is quite small so that 
classical evolution cannot be expected to be reliable after the gluons have been freed and 
one might even worry about the validity of the freeing calculation being done at a classical 
level. 

4.2. Classical fìeld theory versus kinetic theory 

When f g ~ 1/a kinetic theory does not apply and classical field theory is the right 
approach. Thus in the earliest stages of a heavy ion collision, during the time of the 
freeing of the gluons and immediately after, classical field theory is the naturai dynamics to 
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use. There may be a practical problem, however, because the Krasnitz-Nara-Venugopalan 



calculation |22| suggests that there is a rapid lowering of the occupation number as the 
gluons are being liberateci. Thus, although there is a 1/a factor in (28) the actual value 
of f g is not large in any practical circumstance. If it shoulcl turn out that the classical 
calculation is not valici numerically during the freeing of gluons the correct dynamics 
would then be a full quantum field theoretic calculation starting from just before the 
collision up to the time gluons have been freed, a very difficult calculation. 

After gluons have been freed and in the regime (perhaps purely a formai regime) where 
1 < f g < 1/a classical field theory and the Boltzmann equation are equivalent procedures 
for calculating the time evolution of the system ^ . 

When / ~ 1 classical field theory, which only keeps the f g in statistical factors (1 + f g ) 
is an underestimate of the interation strength of gluons while the Boltzmann equation is 
systematic. 



The "bottom-up" picture of evolution ||25|| uses the Boltzmann equation, including elastic 
and inelastic collision terms, and follows the gluon system from shortly after gluons have 
been produced up to the time when equilibration occurs. 



4.3. "Bottom-up" picture []25|];Kharzeev, Levin, Nardi picture ||26[ -p8| 



The "bottom-up" picture of thermalization applies to a collision of heavy ions where 
the freed transverse energy is dominated by gluons having transverse momenta well into 
the perturbative scale. At RHIC energies the transverse momenta which dominate are in 
the region of one GeV and so one might expect the picture to be qualitatively and even 
semiquantitatively correct but precise calculations are always going to be hard to come 
by. The picture does not require that the initial state of the system have large occupation 
numbers but we shall describe it in the context of the initial conditions being given by 
gluon saturated heavy ion-wavefunctions. 

Consider a head-on high energy heavy ion collision. Our description will refer to gluons 
in the centrai unit of rapidity in the center of mass system. The time over which the 
gluons are freed is r ~ 1/Q S . This is also the time scale on which the higher rapidity 
gluons separate, in the longitudinal spatial direction, from the gluons in the centrai unit 
of rapidity with which we are concerned. We presume that a finite fraction, c, of the 
gluons in the initial wavefunction have been freed in the collision as discussed in sec.2.3 
giving initial gluon occupation numbers on the order of 1/a. 

As Q s t grows the occupation number decreases due to elastic scattering during the lon- 
gitudinal expansion of the system. When f g is significantly less than 1/a the Boltzmann 
equation, the basic dynamics in the bottom-up picture, begins to be applicable. When 
Q s r gets as large as ct _3//2 the occupation number for gluons becomes less than one. 

When a~ 3 / 2 < Q s r < « _5//2 the dominant process is soft gluon emission from the hard 
gluons. During this period the density of soft gluons (having pj_ ~ ^fa Q s ) is catching up 
with the density of hard gluons (having p± ~ Q s ), although both densities are decreasing 
in time due to longitudinal expansion. 

When Q s t ~ a~ 5 ^ 2 the density of soft gluons reaches that of hard gluons. At later 
times two important changes happen. (i) The soft gluons stimulate further emission of 
soft gluons so that now the number density of soft gluons grow with time in spite of the 
longitudinal expansion. (ii)The soft gluons equilibrate amongst themselves with the hard 
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Figure 5. 



gluons serving as a source of energy causing the temperature of the soft gluons to grow 
with time as T ~ a 3 Q s T. The rate of flow of energy from hard gluons continues to grow 
until, at Q s t ~ a~ 13 / 5 , the hard gluons disappear, having lost ali their energy, and one is 
left with a completely equilibrated system. At later times the temperature decreases r -1 / 3 
as the equilibrated system continues to expand. The various stages in the equilibration 
process are shwon in Fig.5. Finally, it is important to note that it is probably overly 
optimistic to believe that these separate regimes are sharply distinct at RHIC energies. 

Now let us try to simply characterize charged hadron production at RHIC Take a cen- 
trai collision. Then at mid-rapidity the initial gluon density in the McLerran-Venugopalan 



model is 26 



dN in , 

—4- = c2VR 2 - b 2 pxG (29) 
dyd 2 b 

where c is the fraction of gluons liberated. Eq.(29) is easily obtained by integrating 
(8) over k±, using (9), and inserting the gluon freeing factor c. Now at late times the 
experiments measure the rapidity distribution of charged hadrons. From (29) we get 



< -— — — >= -RcxG(x, <Q 2 >)^ 3.8 (30) 
N par t dy 3 

where N part is the number of participants, the factor 2/3 accounts for the fraction of 
charged particles and R is an inelasticity factor as the partonic system evolves to equi- 
librium. (R also accounts for any, presumably small, change in particle numbers as 
quarks and gluons go to hadrons.) The experimental value[Q of 3.8 corresponds to 
y/s = 200GeV. 
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Figure 6. 



Kharzeev, Levin and Nardi take a very simple picture where R = 1 so that |cxG ~ 3.8. 
The KLN analysis was the first detailed analysis of particle production at RHIC using 
the saturation (color glass condensate) picture and so it has great historical importance. 
I believe it is probably difficult to get cxG as large as 5.7. 

The Krasnitz, Nara, Venugopalan picture^] is also a saturation picture where, how- 
ever, one of the primary achievements was to calculate c and find that c ~ 1/2 at RHIC 
energies. Moderately large values of R and xG are necessary to agree with the data. 

The "bottom-up" picture does not give a precise value for R, and of course there is 
no calculation of c or of xG. However, given values of R, c and xG one does get values 
(estimates) for the time at which equilibration occurs and the temperature at that time. 
For R ~ 3, x ~ 1 and xG ~ 2 one gets T equi i ~ 230MeV and T equ u ~ 3.6 fm, values which 
seem reasonable. 

The KLN, KNV and "bottom-up" pictures are ali based on the same ideas, and differ 
only in the focus of the calculation and in the various details. 

5. Surprises, puzzles, connections 

5.1. Gluon liberation, proton-nucleus vs nucleus-nucleus 

Consider gluon emission in quark- nucleus scattering where the coherence time of the 
gluon is larger than the nuclear radius. The process is best viewed in the rest system 
of the nucleus and is shown in Fig.6 in lowest order. The graph in part (a) of the 
figure corresponds to a gluon in the wavefunction of the quark as the quark reaches the 
nucleus. Graph (b) corresponds to gluon emission after the quark passes the nucleus. The 
calculation is straightforward to perform [[?]]. The square of graph (a) liberates ali gluons 
having k\/Q 2 < 1 while the square of graph (b) gives the "shadow" of graph (a) producing 
gluons with a dk\jk\ spectrum for k\/Q 2 < 1. The interference terms do not contribute 
when k\/Q 2 s < 1 and they cancel the direct terms for k 2 ± /Q 2 < 1. The calculation can be 
generalize to an incident proton on a nucleus with a similar result. Thus one gets 



dN g 2 
dy 

corresponding to a gluon liberation factor of 2. The fact that the gluon liberation factor 
is 2 rather than 1 corresponds to the "shadow" term and has a quantum mechanical 
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interpretation analogous to the result of a = 2 • ttR 2 for the scattering of a particle off a 
hard sphere in non-relativisstic quantum mechanics. 

Now, what is surprising is that the Krasnitz, Nara, Venugopalan calculation gives a 
much smaller number, c ~ 1/2. The lack of a "shadow" terni in heavy ion collisions is 
perhaps naturai. Then we are left with trying to understand why c is near to 1/2 rather 
than 1. In the case of proton-A scattering most of the freed gluons coming from the (a) 
terni in Fig.6 have a momentum, k±, in the proton's wavefunction obeying kj_/Q s <C 1, 
parametrically for large Q s . For A- A collisions the freed gluons in the wavefunction have 
k±/Q s on the order of one so perhaps it is not so surprising that some of the gluons 
having k± on the order of Q s are not freed. This indeed, might be expected if it were not 
for the fact that the freeing process, as described in Sec.4.1 creates so much transverse 
momentum. 

5.2. v 2 (p_l) 

One of the strongest indications of an early equilibration of the quark gluon plasma in 
ion-ion collisions at RHIC is elliiptic flow|30||. Detailed calculations using the equations 



of ideal hydrodynamics suggest that thermalization must occur by r ~ 0.6/m in order 
to get a sufficiently large elliptic flow [pl| ,p2|. On the other hand it appears extremely 
difficult to get such early equilibration from perturbative QCD; the interaction strength 
does not seem strong enough to cause equibibration at such early times[P^-|37[]. One can 
see the difficulty very simply. To match RHIC data the thermalized quark-gluon plasma 
at r ^ 0.6 has a temperature about 330MeV corresponding to a typical quark and gluon 
momentum of about IGeV. Perturbative calculations should be reliable to test whether 
gluons having IGeV can be in equilibrium. The mean free path for a IGeV gluon to have 
a collision with momentum transfer 2GeV is on the order of lOfm at T = 330MeV. Thus 
it is extremely difficult to see how 2 GeV gluons can be thermalized at r ~ O.Qfm. On the 
other hand gluons of 500MeV are stili in the process of being produced at r ~ 0.6/m. 



Indeed, Krasnitz and Venugopalan ||2 1|| find that the dominant energy density only becomes 
produced when r ~ 3/Q s . Thus, while it is not impossible that an effective thermalization 
of gluons in the regime 0.5 — IGeV could occur at r ~ 0.6/m it would be puzzling from 
the QCD dynamics point of view. 

5.3. Saturation and pj_ — broadening 

Experiments on \i— pair production in proton-A reactions observed a broadening in 



the transverse momentum, p±, spectrum of the fi— pairj38]]. On the theoretical side one 
expects the change in the width of the p\ — spectrum to be[|39[ 



APÌ = ^ L P xG (32) 

where L is the length of nuclear material the annihilating antiquark passes through, p is 
the nuclear density, and xG is the nucleon gluon distribution. With xG ~ 1 (32) is in 
reasonable agreement with data. Now recalling (9) and (17) we see that (9) and (32) are 
really the same formula [[40|j. Thus we can get another formula for Q 2 S as 
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Q 2 s (b) 




dAp 2 L 



(33) 



dz 



Now take the experimental number 



dAp 2 ± 



0.02GeV 2 /fm 



(34) 



dz 



for the broadening of the fi— pair spectrum at fixed target energies. This gives 



This number is about a factor of 2 smaller than what one might expect at RHIC, but 



comparison of the fi— pair spectrum of proton-proton and proton-nucleus reactions could 
give a good direct measurement of Q 2 while a comparison of the J j ip production transverse 
momentum spectrum could give a good direct measurement of Q 2 . 
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